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Abstract. We present results from four observations of the ac- 
creting X-ray pulsar Vela X-l with the Rossi X-ray Timing Ex- 
plorer (RXTE) in 1996 February. 

The light curves show strong pulse to pulse variations, 
while the average pulse profiles are quite stable, similar to pre- 
vious results. Below 5keV the pulse profiles display a com- 
plex, 5-peaked structure with a transition to a simple, double 
peak above ~15 keV. 

We analyze phase-averaged, phase-resolved, and on-pulse 
minus off -pulse spectra. The best spectral fits were obtained us- 
ing continuum models with a smooth high-energy turnover. In 
contrast, the commonly used power law with exponential cutoff 
introduced artificial features in the fit residuals. Using a power 
law with a Fermi-Dirac cutoff modified by photoelectric ab- 
sorption and an iron line, the best fit spectra are still unaccept- 
able. We interpret large deviations around —25 and —55 ke V as 
fundamental and second harmonic cyclotron absorption lines. 
If this result holds true, the ratio of the line energies seems to 
be larger than 2. Phase resolved spectra show that the cyclotron 
lines are strongest on the main pulse while they are barely vis- 
ible outside the pulses. 
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1. Introduction 
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Fig. 1. Sketch of the Vela X-l system (center of system). The 
dashes mark the approximate orbital positions when the obser- 
vations were made. The movement and the Roche-deformation 
of HD 7758 1 have been neglected. 



Vela X-l (4U 0900—40) is an eclipsing high mass X-ray binary 
consisting of the 23 Mq B0.5Ib supergiant HD 77581 and a 
neutron star with an orbital period of 8.964 d (van Kerkwijk 
et al., 1995), corresponding to an orbital radius of about 1.7 R* 
(Fig. [!]). The system is at a distance of 2.0 kpc (Sadakane et al., 
1985). Due to the closeness of the neutron star and its com- 
panion, the neutron star (M -1.4 Mq, |Stickland et al., jggg ) cretion disk ([Inoue et al., jggjh. For an in-depth discussion of 



but can also suddenly change to less than 10% of its normal 
value (Inoue et al., 1984). The reason for this variability is still 
unknown. It might be caused by changes in the accretion rate 
due to variations in the stellar wind or the formation of an ac- 



is deeply embedded in the strong stellar wind of HD 77581 
(M - 4 10 -6 Mq /yr, |Nagase et al., 1986[ ). Thus, the system 
is a prime candidate for the study of X-ray production via wind 
accretion. The X-ray luminosity is typically ~ 4 1 36 erg/s, 



the system parameters, see ^an Kerkwijk et al. (1995) 



Vela X-l is a slow X-ray pulsar with a period of about 



283 seconds (Rappaport & McClintock, 1975). Despite signif- 
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icant pulse-to-pulse variatio ns, the pulse profile a veraged over 
many pulses is quite stable ( Staubert et al., 198C ). The profile 
changes from a complex five-peaked structure at energies be- 
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Table 1. Details of these observations of Vela X-l with RXTE. Table 2. Rebinning of the HEXTE-data. The factor is the num- 
The on source time corresponds only to the PCA. ber of instrument channels rebinned into a new channel. 



No. 


date 


JD (center of 


orbital 


on source 


Group 


low 


-high 


factor 






observation) 


phase 


time 


1 


20 


- 39 


4 


1 


1996 February 22 


2450135.66 


0.311-0.324 


6048 s 


2 


40 


- 63 


6 


2 


1996 February 23 


2450137.30 


0.492-0.509 


5696 s 


3 


64 


- 103 


8 


3 


1996 February 25 


2450139.17 


0.701-0.717 


5536 s 


4 


104 


- 153 


10 


4 


1996 February 27 


2450140.79 


0.884-0.895 


5820 s 


5 


154 


- 255 


20 



low ^5keV to a simple double peak ed structure at energies 
above ~15keV ( Raubenheimer, 1990). 



The X-ray spectrum of Vela X-l is usually described by 
a power law with an expon ential cutoff at high energies and 
an iro n Ka line at 6.4 keV ( Nagase et al., 1986 ; White et al., 
1983; |Ohashi et al., 1984|). Below ~3keV, a s oft excess is ob- 
served (^ewis^erjL7T992 ; ^an et al., 1994 ). The spectrum 
is further modified by photoelectric absorption due to a gas 
stre am trailing the neu tron star and due to circumstellar mat- 
ter (Kaper et al., 1994). This results in significantly increased 
absorption at orbital phases ^0.5. In addition, erratic increases 
of tih by a factor of 10 are seen at all orbital phases (Haberl 
& White, 1990). At energies above 20keV, evidence for cy- 
clotron absorption features at ^27 keV and ~54 keV has been 
reported from the High Energy X-Ray Experiment (HEXE) on 
Mir ( [Kendziorra et al., 1992] ), fr om Ginga (|Mihara, 1995[ ), and 
from BeppoSAX observations (Drlandini et al., 1998), but to 
date no statistically compelling result on both lines has been 
obtained. 

In this paper we present the results of the analyses of four 
observations of Vela X-l with RXTE. The observations and 
data reduction are described in Sect. ^. In Sect. |] we discuss 
the light curves and pulse profiles. Section ^ is devoted to the 
spectral analysis, including spectral models, phase averaged 
spectra, phase resolved spectra, and a comparison with other 
instruments. Finally, we discuss the results in Sect. BL 

2. Observations and Data Analysis 

We observed Vela X-l with RXTE four times, covering nearly 
a complete binary orbit from phase 0.3 to 0.9. Each obser- 
vation was between 5 ksec and 6 ksec long (see Table [j] and 
Fig. [j]). The observations were carried out between 1996 Febru- 
ary 22 and 1996 February 27. The exact time and duration of 
each observation are given in Table |l| In our analysis we used 
data from both RXTE pointing instruments, the Proportional 
Counter Array (PCA ) and the High Energy X-ray Timing Ex- 
periment (HEXTE). To avoid scattered photons from the Earth, 
we used only data where Vela X-l was more than 10° above 
the horizon. To extract spectra and light curves, we used the 
standard RXTE analysis software FTOOLS 4.0 and software 
provided by the HEXTE instrument team at UCSD. 

The PCA consists of five co-aligned Xenon proportional 
counter units with a total effective area of about 6500 cm 2 and 
sensitive in the energy range from 2 ke V to ~60 keV (Jahoda 
et al., 1996). Light curves with 250 ms resolution and energy 



spectra in 128 pulse height channels were collected. For spec- 
tral fitting, we used version 2.2.1 of the PCA response matrices 
(Jahoda 1997, priv. comm.). Due to the high statistical signif- 
icance of the data, we applied systematic uncertainties to the 
spectra. These were determined from deviations in a fit to the 
Crab Nebula and pulsar spectrum and had values of 2.5% be- 
low 5 keV, 1 % between 5 and 20 ke V, and 2 % above 20 ke V. 

Background subtraction in the PCA is usually done using 
a background model. This model does not work yet with this 
RXTE epoch 1 data ( Stark, 1997 ), and we were forced to use 
a different approach to background subtraction. Since the PCA 
is not turned off during Earth occults an approximation to the 
average PCA background spectrum can be obtained by accu- 
mulating data measured during the occult. For this, data where 
the source was at least 5° below the horizon was used. Because 
the background depends on geomagnetic latitude, which varies 
throughout the spacecraft orbit, it is obviously not possible to 
generate a time variable background estimate from occultation 
data. Therefore, the light curves presented here are not back- 
ground subtracted. To check the validity of our approach in the 
background subtraction we compared spectra obtained during 
Earth occults from later post gain-shift data with background 
spectra estimated for these data. In general, the agreement be- 
tween the estimated and the measured spectrum was good for 
these data. This suggests that our approach will result in an ad- 
equate background subtraction, especially considering the fact 
that our data are completely source dominated. Finally, due to 
uncertainties in the response matrix above 30keV and since the 
high energy channels might be background dominated, we lim- 
ited the PCA energy range from 3 to 30 keV. 

The HEXTE consists of two clusters of four Nal/Csl- 
phoswich scintillation counters, sensitive from 15 to 250 keV. 
A full description of the instrument is given by Rothschild 
et al. (1998). Our observations were made before the loss 
of energy information from one detector, so we have the 
full HEXTE area. Background subtraction is done by source- 
background rocking of the two clusters which provides a direct 
measurement of the HEXTE background during the observa- 
tion such that no background model is required. This method 
has measured systematic uncertainties of < 1% (Rothschild 
et al., 1998). We used the standard response matrices dated 
1997 March 20 and treated each cluster individually in the 
data analysis. We ignored channels below 20, and to improve 
the statistics of individual energy bins, we rebinned the raw 
(^1 keV wide) energy channels as enumerated in Table 0. 
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2.1. Spectral Models 
2.1.1. Continuum Models 

Despite two decades of work, there still exist no convinc- 
ing theoretical models for the shape of the X-ray spectrum 



in accreting X-ray pulsars (cf. Harding, 1994, and references 



therein). Therefore, we are forced to use empirical spectra in 
the modeling process. 

As in previous observations with Ginga and HEXE (Mi- 
hara, 1995; Kretschmaret al., 1997c), we found it impossible to 



fit the spectra with thermal bremsstrahlung or blackbody spec- 
tra, while we were able to fit our spectra with a power law mod- 
ified by a high energy cutoff. The standard version of this cutoff 
{PLCUT henceforth) is analytically realized as follows (White 
etal., 1983): 



C(E) = E- 




for E < E cut 
for E>E cnt 



(1) 



In this model, the high energy cutoff is switched on at the cutoff 
energy, E cut , and the derivative is discontinuous at this point. 
The source spectrum, however, most probably does not have 
such a sudden turnover. Instead, the importance of the cutoff 
might increase steadily from insignificance at low energies un- 
til it dominates the power law above the cutoff energy. If such 
a spectrum is fit with the PLCUT model, the resulting residu- 
als are likely to mimic a line feature at E cu u and might lead to 
erroneous conclusions (Kretschmar et al., 1997t). 




To demonstrate the effects of applying the PLCUT model 
to spectra with a smooth continuum, we simulated a PCA spec- 
trum using a power law with a smooth exponential turnover, 
including an iron line at 6.4 keV and photoelectric absorption. 
The spectral parameters for the simulation were chosen to be 
appropriate for Vela X-l. We then tried to use the standard PL- 
CUT model to describe this simulated spectrum. The result is 
shown in Fig. ^. The fit features a deep absorption line-like 
feature at about 20keV and some sort of additional cutoff (or 
a second absorption line) above 40 ke V. Note that E cut has al- 
most the same value as the energy of the first line-like feature 
and both are strongly correlated. These two line-like features 
might easily be misinterpreted as a cyclotron absorption line 
and its second harmonic. 

It is obvious, therefore, that the explicit form of the con- 
tinuum is crucial for the detection of absorption (or emission) 
lines like those resulting from interactions between Larmor 
electrons and X-rays at the cyclotron energy. But as the real 
continuum is not known the only thing we can do is to avoid 
those continua from which we know that they introduce artifi- 
cial features in the fit residuals. 

We used the so-called Fer mi-Dirac cutoff (FDCO in the 
following) first mentioned by Ianaka (1986), which has a 
smooth and continuous turnover. 



C{E) = E~ 



(2) 



Energy in keV 

Fig. 2. A simulated spectrum based on a power law with 
a smooth cutoff. This spectrum has been fit with the PL- 
CUT model with E cut «20keV. The ratio is given by (data- 
model)/data. 



Due to historical reasons, the parameters of the PLCUT- and 
the FDCO -model have the same names, but we stress that it is 
not possible to compare them as the continua are different. 

Another spectral model avoiding a sharp turnover is the 
Negative Positive Exponential (NPEX) model introduced by 



Mihara(1995) 



C(E) = (Ai • E- ai + A 2 ■ E +a2 ) x exp 



E 



(3) 



where a\, a 2 > 0. This model can fit very complex continuum 
forms, including dips, nearly flat areas, and bumps strongly de- 
pending on the normalization of the power law with the posi- 
tive photon index. If the second photon index is fixed at a value 
of a 2 = 2.0 this model is a simple analytical approximation 
of a Sunyaev-Titarchuk Comptonization spectrum (Sunyaev & 



Titarchuk, 1980). Mihara (1995) found that this restricted form 



of the model was able to describe a large variety of X-ray pulsar 
spectra. 

We fit our spectra with both the FDCO and the NPEX mod- 
els and found that both provided reasonable descriptions of 



the continua (see Sect. 4.3.1). In comparison, the FDCO spec- 



trum appeared to give more reliable results. With the NPEX 
model the normalization of the second power law was often 
set to practically zero by the fit software, resulting in a stan- 
dard power law with exponential cutoff. We therefore decided 
to use the FDCO model for detailed fits. 
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Fig. 3. PCA light curve of observation 1 . The first gap is due to 
the South Atlantic Anomaly and the second to an Earth occul- 
tation. 



Fig. 4. The beginning of the PCA light curve of observation 2 
with an extended period of time where no pulsations are ob- 
served. 



2.1.2. Absorption Line Models 

Two models are generally employed to describe cyclotron ab- 
sorption line profiles. In both cases, the absorption is described 
by an exponential term of the form e~ T< - E \ where r is the op- 
tical depth of the absorption as a function of energy. This term 



is multiplied by a continuum model (Sect. 2.1.1). The first line 
model is a simple Gaussian profile for the optical depth. The 
second has a Lorentzian profile of the form ( frlerold, 1979 ) 



A(WE/E cyc ) 2 

(E - E cyc y + w 2 



In this case, W is the line width, A the depth, and E cyc the 
line energy. The fundamental line and its harmonics may be in- 
cluded by employing multiple absorption terms, with or with- 
out constraining the relative values of the parameters (e.g. the 
relative line energies or depths). The statistics of our observa- 
tions are insufficient to choose between these two models, and 
all fits reported here employ the Lorentzian profile. The quoted 
widths thus refer to the parameter W. 



3. Light curves and Pulse Profiles 

The average PCA count rate in our four observations is between 
1800cts/s in observation 1 and lOOOcts/s in observations 3, 
corresponding to a normal to low flux level. In addition to this 
normal behavior we observed two special events: 

1 . An increase of the count rate in a flare-like event with 
a maximum PCA count rate of about 7000 cts/s during obser- 
vation 1 (see Fig. |^). Due to an Earth occultation, we did not 
observe the end of the flare: after the occultation, Vela X-l was 
in its normal state again. 

2. An interval of about 550 s duration at the beginning of 
observation 2, where the average total PCA count rate (i.e., 
including the background contribution of ^150 cts/s) was be- 
low 400 cts/s (Fig. |j). During this interval, no pulsations were 
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Fig. 5. Energy averaged PCA pulse profiles for observations 1 
to 4 (panels a to d). The profiles were derived by folding with 
a period of 283.0 sec in 128 phase bins over the entire energy 
range. 



observed. After that, the count rate increased and reached its 
normal value of about 1800 cts/s within one pulse period. This 
behavior has been observed before by Inoue et al. (1984) and 
^apshov et al. (1992) , but its cause is still unknown and further 
investigations are needed. 

We derive an average pulse period of 283.4 ± 1.7 s for all 
four observations, consistent with the value reported by BATSE 
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Fig. 6. Observation 1 pulse profiles in 6 energy bands 
files consist of 128 phase bins. 



for this epoch (BATSE Pulsar Team, 1996). In Fig 
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profiles for the whole PCA energy range are shown 

Despite the fact that these pulse profiles were ot 
different orbital phases, they are quite similar, with 
possibly somewhat affected by the low number of pul; 
(about 20). Energy resolved pulse profiles (Fig. js|) illifetr 
complex behavior found by Raubenheimer (1990) : A 
below ~5 keV, a very complex five peak structure 
The main pulse consists of two peaks with a deef 
tween them, while the secondary pulse consists of 
At higher energies this structure simplifies until the j 
file consists of two pulses with nearly the same intens ity 
~15keV. 

4. Spectral Results 

4.1. Phase Averaged Spectra 

In addition to the FDCO model we allowed for photoelectric 
absorption and an iron fluorescence line at 6.4 keV in 3ur spec- 
tral fits. The iron line is thought to originate in cold, circum- 
stellar material where some X-rays are reprocessed. Vhile we 
found a relatively narrow line in observations 1 and 4 it was 
impossible to obtain an acceptable fit with a narrow line for 
observations 2 and 3. Instead we obtained line widths of about 
1.3 keV (obs. 2) to 1.6 keV (obs. 3; cf. Table |). Similar val- 
ues have been found in some EXOSAT observations ( jottwald 
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Fit with the FDCO model, photoelectric absorption, and 
additive iron line at 6.4 keV to the observation 4. The upper 
hows the data and the model for the PCA and for both 
clusters. The lower panel shows the corresponding ra- 
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995), later ASCA measurements did not find evidence 
small line widths. We caution, however, that our broad 
dths might also be due to the finite energy resolution 
detector in combination with response matrix features, 
the gas stream trailing the neutron star, the photoelectric 
ion increased dramatically between the second and the 
observations. Since the lower energy limit of our spectra 
/ we are not able to observe the soft-excess below 4 ke V 
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reported by [Lewis et al. (1992)| and |Haberl (1994) 
Sti 



1, this description of the data is not acceptable, show- 
lificant deviations near 25 keV and 55 keV. We interpret 
f matures as fundamental and second harmonic cyclotron 
ion lines. Therefore we included a cyclotron absorption 
in our fits. 

ially, we coupled the fundamental and second harmonic 
energies and widths by a factor of 2. Since the statistics 
data were inadequate to constrain both the widths and 
of both lines, we also fixed the width of the fundamen- 



to a value of 5keV (Harding, 1991). While formally 



fits for all observations were found, the second har- 
line had a depth of nearly zero. This contradicts the ap- 
sence of a feature near 55 keV seen in the residuals of 
(see Fig. 7). We concluded that the second harmonic 
coupled to the fundamental line with a factor greater 
10. 
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10 100 
Energy in keV 

Fig. 8. Same as Fig. 7, but with two cyclotron absorption lines. 



Table 3. Parameters of the FDCO model with an Iron line, pho- 
toelectric absorption, and cyclotron absorption lines. The sym- 
bols used are: Nn'. neutral Hydrogen column, in 10 22 cm~ 2 , 
-^cut, Efoid- parameters of the FDCO model (cf. eq. |2|), Fe-u: 
width of the Gaussian iron line in keV (the energy of the line 
was fixed at 6.4 keV), E cyc , Depth 1, and Depth 2: Parameters 
of the cyclotron line, the width of which was fixed to 5 keV, 
Factor: Coupling factor between the fundamental cyclotron line 
energy and the first harmonic, DOF: Degrees of Freedom of the 



fit. 

Obs. 1 Obs. 2 Obs. 3 Obs. 4 





9.2 


+0.5 
-0.7 


4.2 


+ 
-0 


6 
6 


22.3 


+ 
-0 


6 
5 


30.1 


+i 
-i 


1 




a 


1.29 


+0.04 
-0.07 


1.67 


+ 
-0 


03 
04 


1.67 


+0 
-0 


05 
03 


0.94 


+0 
-0 


04 
04 


E cu t 


33.4 


+ 1.4 
-0.5 


39.9 


+1 
-1 


1 

2 


37.7 


+1 

-2 


4 



21.0 


+5 
-7 


1 



-Efold 


8.3 


+ 1.1 
-1.4 


6.9 


+0 
-0 


8 
8 


9.1 


+ 1 
-0 


1 
8 


13.3 


+0 

-0 


4 
5 


Fe-cr 


0.61 


+0.10 
-0.08 


1.29 


+0 
-0 


14 
14 


1.62 


+0 
-0 


17 
32 


0.60 


+0 
-0 


10 
10 


E C y C 


24.1 


+0.7 
-0.5 


22.6 


+1 
-1 


6 
5 


21.6 


+0 
-0 


6 
9 


23.3 


+0 
-0 


4 
3 


Depth 1 


0.18 


+0.05 
-0.01 


0.07 


+0 
-0 


03 
03 


0.16 


+0 
-0 


02 
03 


0.14 


+0 
-0 


02 
02 


Depth 2 


0.32 


+0.10 
-0.79 














1.37 


+0 
-0 


28 
24 


Factor 


2.33 


+ OC 

-0.22 














2.38 


+0 
-0 


10 

10 


Constant 


0.78 


+0.00 
-0.01 


0.71 


+0 
-0 


01 
01 


0.72 


+0 
-0 


01 
01 


0.78 


+0 
-0 


01 
01 



X 2 (DOF) 66 (94) 55 (96) 68 (95) 63 (92) 



Therefore we introduced a variable coupling factor as a new 
free parameter. We used the same parameter to specify the ra- 
tios of line energies and widths (Harding 1997, priv. comm.). 
Using this modified model, we found the fundamental line in all 
four observations, however with variable depth (see Table fy. 
The energy of the fundamental cyclotron line is found between 
22 keV and 25 keV. A second harmonic line is seen in observa- 
tions 1 and 4; however, a second line did not improve the fit for 
observations 2 and 3. Where a line was present, it was coupled 
to the fundamental line by a factor of 2.3 to 2.4 (corresponding 
to a line energy of - 56keV) instead of 2.0. Table | summa- 
rizes the results of the fits to all four observations. To test the 



significance of both lines we used the F-test (Bevington, 1969). 
The results are given in Table [|. While the fundamental line is 
very significant in all four observations, this is not the case for 
the harmonic line. Only for observation 4 is it significant. Note 
that all uncertainties quoted in this paper represent 90% confi- 
dence level for a single parameter. 

4.2. Phase resolved Spectra 
4.2.1. Fine Phase Bins 



In order to study spectral variations with the pulse phase, we 
split the pulse in 10 phase bins as shown in Fig. ^| 

With the reduced statistics of these phase resolved spectra, 
we were unable to constrain the presence of a second harmonic 
line. We therefore restrict this analysis to the continuum and 
fundamental line. 



Table 4. Comparison of x 2 and F-test values of all four obser- 
vations for fits without, with one, and two cyclotron absorption 
lines. Q(F) is the probability that the improvement is acciden- 
tal. 



Obs. 1 



Obs. 2 



Obs. 3 



Obs. 4 



no lines 221.2(97) 

lline 67.1 (95) 

^Q(F) 0.0% 

2 lines 65.6 (93) 

=>Q(F) 34% 



600 r 



69.4 (98) 
55.3 (96) 
0.0018 % 



128.4 (98) 
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63.4 (92) 
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Fig. 9. Definition of the ten phase bins for observation 1 . 
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Fig. 10. Variation of the spectral parameters over the pulse of 
observation number 1. The vertical bars indicate the uncer- 
tainty at the 90% level. 



The behavior with phase is the same for all four observa- 
tions. We therefore discuss only observation 1 (see Fig. |To| ). 
The cutoff energy E cu t, folding energy EWd, cyclotron ab- 
sorption energy E cyc , and the width of the iron line show no 
significant variation over the pulse phase. In contrary, the pho- 
ton index and the depth of the cyclotron line are strongly cor- 
related with the pulse phase. The variation of the photon index 



(Fig. 10) describes a spectral hardening during the pulse peaks. 



This hardening is especially significant in the second pulse. The 
depth of the fundamental line varies also significantly over the 
pulse phase. The cyclotron absorption is strongest in the pulse 
peaks, especially the main pulse, while it is quite weak or ab- 
sent outside the pulses. 

4.2.2. Main Pulse, Secondary Pulse, and Off-Pulse 

To study the second harmonic cyclotron absorption line, it is 
necessary to increase the statistical significance of the bins 




a 1 1 A 1 ill i= 

v A it y >F 



• , I I 1 — I — h 

i b 

+ +++ + ++ 



in 



10 



100 



Energy in keV 



Fig. 11. a and b Fit to the spectrum of the main pulse of ob- 
servation 4. The upper panel shows the folded spectrum and 
model with both cyclotron lines. The two lower panels show 
the residuals for the cases a with both and b without cyclotron 
absorption lines. 



above ^50keV. Therefore, we combined the spectra of the ten 
phase bins into spectra of the main pulse, the secondary pulse, 
and the off-pulse according to Fig. ^. The parameters resulting 
from fits to these spectra are shown in Table |5| The fundamental 
cyclotron absorption line is deepest in the main pulse while it is 
relatively shallow in the secondary pulse and barely detectable 
in the off -pulse. 

For the second harmonic line, the results are less straight- 
forward. It is only found in the main pulse spectrum of obser- 
vations 1 and 4. In both cases it is deep and significant at the 
99.9 % level. A fit to the spectrum of the main pulse of obser- 
vation 4 and the corresponding residuals for models both with 
and without cyclotron absorption lines is shown in Fig. |TT| In 
the case of the other two observations the fit does not improve 
significantly with the inclusion of the harmonic line. This is 
almost certainly due to the low statistical quality of the data 
above ^60keV. The harmonic line cannot be observed in the 
secondary pulse and the off-pulse except for observation 4, 
where it is seen in both. 



4.2.3. Pulse minus off-pulse 



The results from Sect. 4.2.2 indicate that the cyclotron absorp- 



tion lines originate primarily in the pulsed emission and not 
in the persistent flux. To further study this effect, we created 
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Table 5. The fit-results of the fits to spectra of the main, the secondary, and the off -pulse. The width of the fundamental cyclotron 
line was fixed to 5 keV and the energy of the iron line to 6.4 keV. The symbols have the same meaning as in Table (§ 



Observation 1 



Observation 2 



Parameter 


Main Pulse 


Sec. Pulse 


Off-Pulse 


Main Pulse 


Sec. Pulse 


Off-Pulse 


iV H [ 10^ cm" 2 ] 


°-- > -0.9 


O- 1 -0.8 


li-f- -0.9 


4 5 +u ' v 

-0.6 


q 8 + u - 4 
->-° -0.5 


3.4 


+0.5 
-0.4 


a 




1 20 + 007 
i - zu -o.os 




i-6i ±S:SI 


i cq +0.05 
-0.05 


1 JTJT +0.02 

-0.03 


1.81 


+0.02 
-0.02 


E C ut 


[keV] 


32.4 ±l 3 6 


33.6 ±S;g 


32.4 I 2 ; 4 


35 1 + 20 

-2.6 


40.7 +^ 


44.9 


+ 1.4 
-1.2 


-Efoid 


[keV] 


11-3 ±ii 


7 6+ 14 

'•° -1.3 


8-0 


74 +0-9 

//f -0.8 


6 +o:s 


4.6 


+ 1.1 
-1.0 


Depth 1 




0-28 ±™ 


n , o +0.04 
U.18 _ .04 


10 + 006 

"■ lu -0.06 


31 + 006 
UJ1 -0.06 










[keV] 


24 3 + - 4 

z * f,J -0.4 


24.9 ±{° 


23.4 tH 


oa a +0.6 
z " , "> -0.5 








Depth 2 




1 55 + - 72 

l.JJ -0.64 














Factor 




2 40 +0 - 21 

-0.13 














Fe-cr 


[keV] 


54 + - 20 

U - J ^ -0.23 


0-58 ±C;1| 


79 + 16 


1 42 + ' 15 

-0.15 


1 ni +0.12 
1J1 -0.12 


1.03 


+0.17 

-0.16 


Const. 




n 7R +0.01 
U./8 _ 002 


0.79 ig;gj 


77 +0.02 

' -0.01 


72 + - 02 

u - IA -0.02 


7S + 001 

' J -0.02 


0.68 


+0.02 
-0.01 



63.9 (94) 76.9 (96) 50.3 (96) 59.0 (96) 



(DOF) 



Observation 3 
Sec. Pulse 



Observation 4 



Parameter 



Main Pulse 



IVh [10^ 2 cm" 2 ] 
a 

Scut [keV] 
E ioM [keV] 
Depth 1 

E cyc [ keV] 
Depth 2 
Factor 

Fe~cr [ keV] 
Const 



+ l.b 

-1.2 

+0.11 

-0.15 

+7.0 

-14.0 

+2.6 

-2.3 

+0.06 

-0.06 

+0.7 

-0.7 



?! 3 + 11 
zl -' -1.0 

1 AT, +°- os 
i-'+J -0.09 

32.9 ±H 
10.1 ±\i 

o.io ±° ° i 

20.9 ±\ 6 6 



Off-Pulse Main Pulse Sec. Pulse Off-Pulse 



"+m3 — 

-0.9 

+0.06 

-0.06 

+ 1.9 

-2.4 

+2.6 

-2.0 

+0.04 

-0.03 

+ 1.5 

-1.4 



+ 1.0 
-1.1 

+0.03 

-0.07 

+ 1.7 

-0.8 

+0.7 

-0.6 

+0.03 

-0.04 

+0.4 

-0.3 

+0.67 

-0.61 

+0.15 

-0.15 

+0.50 

-0.48 

+0.01 

-0.02 



+ 1.2 

-1.0 

+0.04 

-0.05 

+ 1.7 

-0.8 

+ 1.0 

-0.9 

+0.03 

-0.03 

+ 1.0 

-0.9 

+0.34 

-0.35 

+0.17 

-0.11 

+0.16 

-0.16 

+0.01 

-0.02 



+ 1.5 

-1.7 

+0.10 

-0.13 

+5.9 

-9.2 

+2.0 

-2.0 

+0.03 

-0.04 

+0.7 

-0.6 

+0.75 

-0.36 

+0.20 

-0.12 

+0.11 

-0.12 

+0.01 

-0.01 



22.5 
1.57 
25.9 
12.2 
0.20 
21.5 



1.51 
0.74 



1.61 



+0.23 
-0.23 

Hg-.gl 0-72 
79.8 (96) 



21.6 
1.85 
44.7 
6.4 
0.11 
19.9 



+0.20 , aq +0.17 
-0.19 l Hy -0.17 



+0.02 n 7n +0.02 
— 0.01 u - /u — 03 



35.8 

I. 28 
31.5 

II. 2 
0.29 
23.7 
1.41 
2.52 
0.80 
0.77 



29.2 
0.76 
27.2 
9.6 
0.16 
24.2 
0.76 
2.12 
0.46 
0.79 



30.5 
1.27 
29.7 
13.9 
0.10 
23.4 
1.49 
2.35 
0.59 
0.78 



43.3 (96) 63.5 (96) 76 (94) 60.5 (94) 57.9 (94) 



X" 



(DOF) 



Table 6. x 2 an d F-test values for fits to spectra of the main, the secondary, and the off-pulse without, with one, and with two 
cyclotron absorption lines. 



Observation 1 



Observation 2 





Main- 


Sec- 


Off- 


Main- 


Sec- 


Off- 




number of lines 


DOF 


Pulse 


Pulse 


Pulse 


Pulse 


Pulse 


Pulse 





(98 DOF) 


302.3 


132.4 


58.8 


140.6 


87.6 


78.9 


1 


(96 DOF) 


74.7 


76.9 


50.3 


59.0 






<2(F) 




0.0% 


< 10" 9 % 


0.06 % 


0.0% 






2 


(94 DOF) 


63.9 












Q(F) 




0.06 % 



















Observation 3 






Observation 4 








Main- 


Sec- 


Off- 


Main- 


Sec- 


Off- 


Number of lines 


DOF 


Pulse 


Pulse 


Pulse 


Pulse 


Pulse 


Pulse 





(98 DOF) 


131.2 


68.2 


93.1 


339.7 


152.2 


100.2 


1 


(96 DOF) 


79.8 


53.2 


63.5 


87.6 


73.7 


62.5 


<2(F) 




10" 8 % 


0.0003 % 


< 10" 5 % 


0.0% 


10" 8 % 


10" 8 % 


2 


(94 DOF) 








78.7 


60.5 


57.9 


Q(F) 










0.065 % 


0.1% 


2.8% 



a "pulsed" spectrum by accumulating the main or secondary 
pulse spectra and subtracting the off-pulse spectra. The results 
of fits to these spectra are shown in Table 0: the fundamental as 
well as the second harmonic line (where observable) are deep- 



est on the main pulse and less deep, but still significant, in the 
secondary pulse. 

The energy of the cyclotron absorption lines do not vary 
significantly between observations. From the pulsed spectra we 
derive average line energies of 24.2 keV for the main pulse 
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Table 7. Parameters of the cyclotron absorption lines in the 
"pulsed" spectra for all four observations. 



Table 8. Comparison of the NPEX and Lorentzian cyclotron 
line parameters for phase average spectra found by Mihara 



Main Pulse 


Param. 


Obs. 1 


Obs. 2 


Obs. 3 


Obs. 4 


Depth 1 
Depth 2 
Factor 


69+ 06 

7 91 +1.07 
z - zi -1.63 

Z - JJ -0.12 


24 6 +L ' 2 

z ^-°-0.8 

1 04+ - 57 


23.8 + i* 
3<j+ - 12 

U - JJ -0.12 

57+ - 83 

u -"-0.57 

2.15i~ 


24 +u;/ 

z ^- u -0.6 

65+° 12 
u - DJ -o.io 

1 63+ 1 ' 33 

9 «+0.69 
i - JJ -0.28 






Secondary Pulse 




Param. 


Obs. 1 


Obs. 2 


Obs. 3 


Obs. 4 


Ecyc 

Depth 1 
Depth 2 
Factor 


70 ,+o.y 

^J-J-0.9 

ua '-0.09 
, 7 q+0.65 

24 n+0.19 
z -^~-0.13 






24.7+i' 3 

o n+ 04 

u - 1:, -0.04 
f) Q4.+0-35 
U - y f--0.35 

2 10+° 15 



and 24.0 keV for the secondary pulse, the difference being well 
within the uncertainty. This seems to be even more the case as 
there does not seem to be any correlation between the pulse 



phase and the line energy in the ten phase bins (see Sect. X.2. 1 
and Fig. [[(]). 

The coupling factor between the fundamental and second 
harmonic line is greater than 2 in all four observations for those 
pulse phases where a second line is observed. We derive an 
average coupling factor of 2.4 which agrees with the values 
previously obtained from the phase averaged spectra and the 
main pulse spectra. 

4.3. Comparison with other instruments 
4.3.1. Ginga 

The Ginga -team used the NPEX model modified by photoelec- 
tric absorption and an additive iron line at 6.4 keV (Mihara, 
1995). They found a cyclotron absorption line at 24keV and its 
second harmonic line. They used a fixed coupling constant of 
2.0. 

In order to compare our results with the Ginga observa- 
tions, we performed fits using the NP EX model. Table [| shows 
the values found by Mihara (1995) using Ginga data and fits 
with the NPEX model to our data. Finally, we compared this 
fit to fits with the FDCO model to see which model describes 
the data better. 

Table [| shows that our results are in very close agreement 
with the Ginga observation of the fundamental cyclotron ab- 
sorption line. Both the FDCO and NPEX models provide very 
similar line parameters. The only significant difference is that 
we find a coupling factor of about 2.3 (independent of the 
model), while Mihara used the canonical fixed value of 2.0. 



4.3.2. HEXE and TTM onboard Mir 

By using both instruments in their work, Kretschmar et al. 
(1997c) were able to analyze the broad energy range from 2 
to 200 keV, i.e., very similar to RXTE but lacking the effective 
area especially in the lower energy band. They used the stan- 



(1995) with Ginga, Orlandini et al. (1998) with BeppoSAX 



and our observation 4 RXTE data. We also show a fit to our 
data with the FDCO model. For the Ginga and RXTE data, the 
width of the second harmonic is related to that of the funda- 
mental by "Factor." 



Param. 



Ginga 



BeppoSAX 



RXTE 



NPEX 



FDCO 





u ' ol -0.05 


34+o.n 

U-Jt_o.ll 


42 +o.02 

u -^ z -0.02 


94+U-U4 

u -^-0.04 


OL1 


-2 fixed 


n I +0.5 
z - 1 -0.5 


-2 00+ ' 04 

z - uu -0.08 




-Efold 


6 4+0-1 
°-^-o.i 




7 o+0.04 

'•-'-o.os 




Ecyc 


24 <> +0 ' 5 

z ^- 3 -0.5 


CQ+2.0 
-"-1.0 


23 4 +0 - 4 

z ->-^-0.3 


70 O+0.4 
z - , --'-0.3 


Depth 1 


u.uoj_ 015 


1 5+ ' 7 

1 --'-0.6 


ul '-0.02 


14+ 002 

ua ^-0.02 


Width 1 


2 2+ 1 - 


201* 


5 fixed 


5 fixed 


Depth 2 


o.8ot°;| 




1 00+ 008 

1 - uu -0.07 


, 07+0.28 

l-->'_0.24 


Factor 


2 fixed 




7 oo+0.09 

^■"-o.os 


7 oQ+0.10 



dard model: a power law with high energy cutoff, photoelectric 
absorption and an iron line. They found a cyclotron absorption 
line at 22.6 keV very similar to our value. They also find a deep 
second harmonic absorption line which is coupled to the fun- 
damental line by a the fixed factor of 2.0. Fits to the HEXE 
and TTM data using the FDCO model (Kretschmar 1997, priv. 
comm.) resulted in photon indices, cutoff energies, and folding 
energies which are very similar to ours. 

4.3.3. BeppoSAX 

Orlandini et al. (1998) analyzed Vela X- 1 spectra obtained with 
BeppoSAX. They simultaneously fit multi-instrument, broad- 
band (2-100keV), phase-averaged spectra using the NPEX 
and Lorentzian cyclotron absorption models. The resulting pa- 
rameters are given in Table ^|. While the BeppoSAX NPEX 
parameters agree well with both Ginga and RXTE, no line 
near 25 keV was required to achieve an acceptable fit. Based 
on this fact and based on the BeppoSAX Phoswich Detector 
System (PSD) data from 20 to 100 keV alone, they argued that 
the line at ^54keV is the fundamental one. However, they re- 
port that allowing 2 absorption lines in their fit, the first con- 
strained to lie between 10 and 40keV, results in a weak feature 
at ^24keV. While this result was not statistically significant, 
it makes clear that their data can admit a feature near 25 keV. 
Further, their broad-band fits are based on phase-averaged data 
only, while the present work and phase-resolved Ginga (Mi- 
hara, 1995) analyses show that the line is strongest during the 
pulse peaks. According to the results in Table |[ the ^25keV 
line is quite weak in the phase average, having an optical depth 
only near 0.1. Given that HEXE /TTM, Ginga, and now RXTE 
all detect a line near 25 keV, it seems most likely that this is the 
fundamental energy, corresponding to a surface field of about 
2 10 12 G. As prlandini et al. (1998)] note, it is also possible that 
the line strength is variable with time, making the lack of a 
detection with BeppoSAX less unlikely. 
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5. Discussion and Summary 



the 



The RXTE observations of Vela X-l demonstrate 
ity of this source on many time scales. Besides the 
strong pulse to pulse variations and the varying intehsity 
we observed a flare and an interval with no obsen 
pulsations. 

In contrast, the fundamental accretion and emiss 
try is obviously very stable, as the comparison of 
pulse profiles with i n our observations and with 



variabil- 
i|vell known 
states 
able X-ray 



'tie 



o her 



flWhite et al„ 1983fc [Nagase et al„ 198$ |Choi et al 



landini et al., 1998) demonstrates. Our energy 
profiles show the transition from the complex five 
ture at lower energies to the double peak above 15 
precedented detail. The determined pulse period 
tinuation of the current overall spin-down trend. 

One of the main goals of our observations 
or disprove the existence of cyclotron absorption 
around ^25 ke V and ~50 ke V has turned out to be 
cult than anticipated. This is due to the fact that th 
the conti nuum model spectrum affects the results si 
As Sect. O 



resolved pulse 
3eak struc- 
keV in un- 
a con- 



shaws 



li le 



demon strates, the unwary use of the 
high energy cutoff ( White et al., 1983 ) with its 
of the cutoff can introduce line-like residuals and 



the introduction of artificial features in the modeling 



to confirm 
features 
more diffi- 
; choice of 
nificantly. 
'classical" 
onset 
thus lead to 
This fact, 



al irupt 



al., 1997a 



which we have reported previously ( Kretschmar et 
Kretschmar et al., 1997b), has led us to study mjodels 



a smooth turnover, especially the FDCO and 
(Eqs. ||and |]). Both fit the data equally well, so w( 
centrated on the FDCO model which has one 
and somewhat less freedom in its overall shape 

Still - regardless of the continuum description 
could only achieve a convincing fit to the phase av 
tra by including a fundamental cyclotron line at 22 
our models. According to the F-test (Table |j) thif 
significant for all observations. We therefore conch de 
fundamental line energy is 20-25 keV in contrast to 
of Drlandini et al. (1998). For two of the four obs 



used - we 
aged spec- 
:!4keVinto 
feature is 
that the 
the results 
;rvations a 

harmonic feature is also observed, but with lesser significance. 
Contrary to the usual assumption, the second harmonic line 
seems not be coupled by a factor of 2.0 in centroid energy and 
width to the fundamental, but by a factor of 2.3 to 2.5 instead. 
Note, however, that we can not formally rule out a factor of 
2.0 and a weak second harmonic due to limited fit Statistics at 
higher energies (cf. Sect 



Confirming earlier results (Kretschmar et al., 



hara, 1995), th e pulse phase resol ved spectra indicate that, as 
with Her X-l ( |Soong et al., 1990j ), the line features are deep- 
est and most significant in the main pulse and barel) detectable 
outside the pulses. Again, the best fit results are ob ained with 
a coupling factor of >2 between fundamental and i econd har- 
monic. 

At the moment we have no firm explanation 
prising result. Relativistic corrections will cause the 
energies of the lines to decrease, with increasing po 
the incident photons relative to the magnetic field 



on geome- 
averaged 
results 



1996: Or 



with 

NP$X models 
have con- 
partmeter less 



997c; Mi 



f >r 



this sur- 
resonance 
ar angle of 
( |Harding & 



aid 



Daugherty, 1991) 
ing harmonic nu nber 
< 2. The possit ility 
ferent regions at 
fact that we obs^r 
resolved s 

If the emission 
the "fan-beam' 
could in principl 
ferent effective 
is not clear how 
pulse structure 
sion geometry at 

Another poss ibility. 
our results are 
models and line 
Harding, 1996) 
clotron lines dev 
used throughout 
rently no more 
tral analysis. 

Although we 
presented here 
tent models to 
ments deliver. 



region is significantly extended in height - 
cenario - the main contribution to the lines 
e be at different heights and therefore at dif- 
strengths (cf. |Burnard et al., 199 1| ). But it 
this can be reconciled with the simple double 
therefore probably relatively simple emis- 
these energies. 

', which needs further exploration, is that 
distorted by assuming too simple continuum 
shapes. Monte-Carlo calculations (Araya & 
lemonstrate that the expected shapes of cy- 
ate strongly from the simple Lorentzian shape 
in spectral analysis. But, unfortunately, cur- 
cibmplex model exists in suitable form for spec- 



t ler 

fully 
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